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We describe a fast direct numerical simulation (DNS) method that promises to di¬ 
rectly characterise the hydraulic roughness of any given rough surface, from the hy¬ 
draulically smooth to the fully rough regime. The method circumvents the unfavourable 
computational cost associa ted with sim ulating high-Reynolds-number flows by employ¬ 
ing minimal-span channels ( Jimenez fc Moinlll99P . Proof-of-concept simulations demon¬ 
strate that flows in minimal-span channels are sufficient for capturing the downward 
velocity shift, that is, the Hama roughness function, predicted by flows in full-span chan¬ 
nels. We consider two sets of simulations, first with modelled roughness imposed by body 
forces, and second with explicit roughness described by roughness-conforming grids. Ow¬ 
ing to the minimal cost, we are able to conduct DNSs with increasing roughness Reynolds 
numbers while maintaining a fixed blockage ratio, as is typical in full-scale applications. 
The present method promises a practical, fast and accurate tool for characterising hy¬ 
draulic resistance directly from profilometry data of rough surfaces. 


Key words: Authors should not enter keywords on the manuscript, as these must be cho¬ 
sen by the author during the online submission process and will then be added during the 
typesetting process (see http://journals.cambridge.org/data/relatedlink/jfm-keywords.pdf 
for the full list) 


1. Introduction 

Scientists have, for years, been documenting the relationship between surface rough¬ 
ness and its hydraulic resistance, the former pertaining to geometry of the surf ace w hile 
the latter is dictated by the dynamics of flow over the roughness (see reviews bv ljimene3 
20041 : Flack fc Schultz) 201 0). The cataloguing process can be a never-ending exercise be¬ 
cause the characteristics of each rough surface are unique. In order to make predictions 
on full-scale operations, it is necessary to establish the equivalent sand-grain roughness 
k s of a given surface, which relates the drag increment of the given surface to an equiv¬ 
alent surface composed of monodisperse sand grains. Once k s has been determined, it 
is possible to pred ict the drag penalty at application Reynolds numbers using either 
the Moody chart ( Moodvl [l944h . for pipe and channel flows, or developments of thi s 
for zero-pressure-gradient boundary layers ( Prandtl &; Schlichtind 19551 : Granvillelll958 ). 
Generally speaking, the approach to date has been to first identify a particular rough sur¬ 
face of scientific or engineering interest, and then to characterise its hydraulic resistance 
through well-controlled laboratory experiments. By exposing the rough surface to an in¬ 
creasing range of flow speeds, until such point beyond which the resistance coefficient Cf 
becomes constant, referred to as the ‘fully rough’ asymptote, it is possible to ascertain 
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k s . A convenient dimensionless group here is the equivalent roughness Reynolds number 
= k s U T /v , where the friction velocity U T = yV 0 / p = U x ^/Cj/2] r 0 is the wall drag 
force per plan area; p is the density; Uoo is the freestream velocity or centreline velocity 
for internal geometries; and v is the kinematic viscosity. For accurate predictions, it is not 
sufficient to merely establish k s , rather the increment of Cf caused by surface roughness 
must be mapped with respect to at all conditions from the dynamically smooth up 
to the fully rough regimes, which generally covers the approximate range 5 < k + < 100. 
Over the last century or so, this painstaking and time-consuming procedure has been 
repeated for many surfaces of interest and as a result, there is now a rich database of 
roughness amassed in the published literature over time. 

One problem here is that there is no widely applicable function that relates k s , a 
dynamic parameter, to readily observable or measurable geometric properties of the 
surface, such as the root-mean-square roughness height k rms or average roughness height 
k a . Certain classes of surfaces, say sand-grain roughness, may exhibit an approximate 
proportionality between k s and some physical surface length scale, but as a general rule 
this proportionality will not hold across different surface topologies. There are numerous 
attempts in the literature to formulate more complicated functions, often involving some 
measure of mean surface height such as k rms or k a and other properties relating to the 
shape or arrangeme nt of roughness elemen ts such as skewness, effective slope, solidity and 
so on (see review bv lFlack fc Schult3l201(j) . Though many of these parameters have some 
success in describing the particular class of surfaces for which they were formulated (e.g. 
painted or sanded surfaces), none are widely applicable across the almost limitless range 
of surface topologies that are encountered in engineering and meteorological applications. 

The present method to directly evaluate k s represents a paradigm shift from such 
parameterisations that are based on a handful of geometrical factors. Recognising that all 
roughness geometries are unique, and that a one-size-fits-all formulaic solution is proving 
elusive, we have sought an approach to minimise the expense involved in experimentally 
determining k s . Specifically, the present method relates raw profilometry data, which 
is described by many degrees of freedom, directly to k s . In contrast, previous indirect 
approaches first reduce raw profilometry data to a limited subset of geometrical factors 
such as k rrns and skewness before relating these factors to k s . The inevitable loss of 
information in reducing raw profilometry data to geometrical factors is avoided by the 
present method and, in this sense, the present method represents a direct evaluation of 
the roughness function. 

The approach relies on the direct computation of hydraulic resistance by direct nu¬ 
merical simulation (DNS), which we presently show can be made substantially cheaper, 
and therefore faster, than previously thought. We circumvent th e otherwise-pro hib itive 
cost of conve ntional DNS by employing minimal-span channels ( Jimenez fc MoiiilllQOll : 
Hwang 2013), the rationale of which is discussed in the following. Normally, a straightfor¬ 
ward and direct computation of roughness drag using DNS employing full-span channels 
is extremely expensive as it entails simultaneously capturing both the bulk flow, which 
scales with the half-channel height, h , and the near-wall flow around the roughness ele¬ 
ments, which scales with a chara cteristic roughness height, k. Given a rough s urface of 
fixed blockage ratio k s /h < 1/40 ( Jimened 2004 : Flack. Schultz fc Shapirol[2005h . a com¬ 
plete characterisation of hydraulic resistance requires parametric simulations that sweep 
through the equivalent roughness Reynolds numbers, kf ~ 5 to 100, corresponding to 
the hydraulically smooth and the fully rough regimes, respectively. For the blockage ra¬ 
tio, h/k s = 40, this means performing parametric simulations at the friction Reynolds 
numbers, Re T = hU T /v = (h/k s )(k s U T /v) « 200 to 4000, which are currently unfeasible. 
Recall that the cost of DNS, counting the number of spatial and temporal degrees of free- 
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dom, scales unfavourably as Re 3 (Ipondliooil §9.1.2). For grids conforming to the surface 
of the roughness elements, this cost is further exacerbated by the need for increased mesh 
density, and reduced time steps. 

However, the extreme cost associated with conventional DNS employing full-span chan¬ 
nels seems unnecessary. The quantity of interest from an engineering point of view is the 
retardation in the mean flow over the roughness relative to the smooth-wall flow. This 
relative flow retardation or downward velocity shift, A U, occurs mostly in the vicinity of 
the roughness layer, but holds constant above a few roughness heights, well into the log 
layer (if it exists) and the wake regi on, a behaviour described by Townsend’s outer-layer 
similarity hypothesis ( Townsend 1976ll . This suggests that a simulation of only the near¬ 
wall region and its interaction with the roughness geometry is required in order to extract 
A U + = AU/Ur , which is known as the (Hama) roughness function. This distillation of 
the problem is consistent with the observation that A U + does not depend on the bulk 
flow but only on kf and other details of the roughness geometry. 

A framework for simulatin g onl y the near -wall dynamics is the minimal channel. The 
concept is first described bv ljimenez fe Moinl ( 199lh and is currently rec eiving renewed at¬ 
tention in various c o ntexts of understanding wall-bou nded turbulence ( Flores fe Jimene3 
2010t Hwand 20131 Lozano-Duran fe Jimened 20141) . Presently, we exploit this frame¬ 
work for measuring AU + by fully resolving the near-wall Navier-Stokes dynamics and 
its interaction with the roughness geometry. The prohibitive cost of conventional DNS 
is alleviated by use of these minimal-span channels, which are designed to preclude the 
bulk flow that scales with h. Without the bulk flow, the cost of DNS with roughness 
now only scales as fc+ 3 , which is quite feasible for the engineering task at kf ~ 5 to 
100. In principle, the computational cost is potentially ( h/k s ) 3 times less than that of a 
conventional DNS in a full-span channel. 

In the remainder, we demonstrate the efficacy of this appro ach and develop gu i deline s 
for its use. Beginning first with the parametric forcing model of Busse fc Sandhaml ( 2012 ^. 
we carefully confirm that the minimal channels return the same estimate for the roughness 
function as that given by full-span simulations (§ EH). An important subtlety here lies 
in the choice of the sminimal spanwise unit, which we demonstrate is related not only 
to the usual span of the near-wall cycle (approximately 100 viscous wall units), but 
also to some physical roughness height k and, presumably, the roughness texture under 
investigation. Having established the accuracy with the simple forcing model, we then 
test this procedure with the more realistic test case of a three-dimensional roughness, 
explicitly described by grids conforming to the rough surface, and compare results with 
full-domain simulations (§& 


2. Direct numerical simulations 

We present two sets of DNSs (table [lj). In Set 1, roughness is mo delled by body 
forces (figure [T|c) in a fourth-order staggered-grid code as described bv lMorinishi et al. 
( 19981 ). The code is wr itten by the first author and has been used in other studies 
Chung et ali 12014 ). In Set 2, roughness is explicitly r epresented by a roug hness- 


l e -g- - , . - _ .. 

conforming grid (figure EK) in the code CDP described bv iMahesh et all ( 20041 ). The 
purpose of the Set-1 simulations with modelled roughness is to test the present minimal- 
channel method for a hypothetical roughness in the absence of any imposed wall-parallel 
length scales. Such a configuration, which we shall call homogeneous roughness, is ideal 
for investigating how the characteristic roughness height (as opposed to wall-parallel 
length scales) sets the minimum allowable span L y of the channel such that AU + is 
still accurately captured. The purpose of the Set-2 simulations with grids conforming to 
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Figure 1 . Computational domain of (a) the full channel and (6) the minimal channel. The 
magnified view of the channel surface illustrates the rough wall represented by (c) the rough¬ 
ness-forcing region, 0 < a < k, or (d) the explicit sinusoidal roughness with semi-amplitude k 
and mean a,t z — 0. Both top and bottom channel walls are rough. 


the rough surface is to demonstrate that the present minimal-channel method works as 
expected for a specific roughness geometry, namely th e sinusoida l ‘ egg-c arton’ roughness 
that has previously been characterised (in a pipe) bv IChan et al. ( 2015 1. 

In both sets of simulations, we solve the following Navier-Stokes equations of motion: 


dut d{ujUi) 
dt dxn 


1 dp 
p dx 


d 2 


du i 




dx] 


dxj 


(2.1a,6) 


where it,; is the velocity; t is time; Xj is the spatial coordinate; p is the pressure; and 
f(t) is the spatially uniform, time-varying, mean pressure gradient that drives the flow 
at constant mass flux. The streamwise, spanwise and wall-normal directions are referred 
to as either x±, X 2 and X 3 or x, y and z respectively. Periodic boundary conditions are 
imposed in the streamwise and spanwise directions with the respective domain sizes, L x 
and L y (figure [Da, b). In the reference full-span simulations, L y measures irh for Set 1 
and 27 rh for Set 2 (at several selected Re T ), while in the minimal-span simulations, L y is 
constrained by the consideration of three issues: 

(1) the span of the near-wall cycle (assumed to be approximately 100 v/U T ), 

(2) the span of the roughness elements, \ y (absent for Set-1 simulations), and 

(3) the physical height of the rough ness elements, k. _ __ 

Regarding point (1), previous studies ( Jimenez fc MoiiJ 199lfc iHwan 3HoH) have shown 
that L+ > 100 or else the near-wall cycle, which is the signature of wall turbulence, 
cannot be properly captured. Regarding point (2), the span of the channel must be wide 
enough to accommodate the widest scale of the roughness elements in order to properly 
represent the nature of the roughness in question. For the Set-1 simulations, which have 
a homogeneous forcing, this point is not an issue. Regarding point (3), it turns out that 
the height of the roughness elements also plays a role in setting the minimum L y , an 
issue that will be discussed in § 13.II T he chosen reso l utions (table [D ) are co mparable 
to other channel-flow DNSs (e.g. 


Moser et all 1999t Bernardini et al. 2014 1 and the 


streamwise domain lengths are long enough to accommodate several near-wall streaks, 
which are approximately 1000 1 '/U T . Although the present study focusses on minimising 
the spanwise dimension of the channel, w hich is the critical dimension in setting the 
dynamics of the log layer ( Flores fc Jimene~3l2010ll . one presu mes that the strea mwise 
dimension will also be subject to constraints. For a smooth wall, Chin et all ( 2010h show 
that the mean velocity profile is slightly overestimated when L+ < 1000 (at Re T ss 180 for 
pipe flow), suggesting that the near-wall streaks need to be properly captured. It is also 
clear that L x should be greater than X x , where is the largest characteristic streamwise 
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Roughness 

Span 

Re T 

h/k 

L+ 

L t 

N x 

Ny 

N z 

Ax + 

Ay+ 

A z+ 

A*+ 

’ Modelled 

Minimal 

180 

40 

3707 

116 

384 

24 

192 

9.7 

4.8 

0.02 

2.9 

Modelled 

Minimal 

395 

40 

3707 

116 

384 

24 

192 

9.7 

4.8 

0.05 

6.5 

Modelled 

Full 

590 

40 

3707 

1854 

384 

384 

256 

9.7 

4.8 

0.04 

7.2 

Modelled 

Minimal 

590 

40 

3707 

116 

384 

24 

256 

9.7 

4.8 

0.04 

7.2 

Modelled 

Full 

950 

40 

5969 

2985 

640 

640 

384 

9.3 

4.7 

0.03 

7.8 

Modelled 

Minimal 

950 

40 

3581 

112 

384 

24 

384 

9.3 

4.7 

0.03 

7.8 

Modelled 

Minimal 

2000 

40 

3707 

116 

384 

24 

768 

9.7 

4.8 

0.02 

8.2 

Modelled 

Minimal 

2000 

40 

3707 

232 

384 

48 

768 

9.7 

4.8 

0.02 

8.2 

Modelled 

Minimal 

2000 

40 

3707 

463 

384 

96 

768 

9.7 

4.8 

0.02 

8.2 

Modelled 

Minimal 

4000 

40 

3707 

116 

384 

24 

1024 

9.7 

4.8 

0.02 

12.3 

Modelled 

Minimal 

4000 

40 

3707 

232 

384 

48 

1024 

9.7 

4.8 

0.02 

12.3 

■ Modelled 

Minimal 

4000 

40 

3707 

463 

384 

96 

1024 

9.7 

4.8 

0.02 

12.3 

Explicit 

Full 

180 

18 

2262 

1131 

510 

254 

156 

4.4 

4.5 

0.13 

4.5 

Explicit 

Minimal 

180 

18 

4029 

141 

910 

32 

156 

4.4 

4.4 

0.13 

4.4 

Explicit 

Minimal 

360 

18 

4100 

141 

924 

28 

314 

4.5 

5.0 

0.13 

5.0 

Explicit 

Minimal 

540 

18 

4029 

212 

755 

42 

470 

5.2 

5.0 

0.13 

5.2 

Explicit 

Minimal 

1080 

18 

4241 

424 

869 

84 

940 

4.8 

5.0 

0.13 

5.0 


Table 1. Simulation cases listed with nominal Re T . Each of these 17 cases is run using both 
smooth and rough walls, that is, a total of 34 separate simulations are run. In Set 1 (modelled 
roughness) ak = 1/40 and in Set 2 (explicit sinusoidal roughness) \ x /k = \ y /k « 7.1. The 
grid spacing is uniform in the streamwise and spanwise direction, while the grid spacing in the 
wall-normal direction is stretched, with the finer wall grid spacing Az w and the coarser centreline 
grid spacing A z c . 


scale of the roughness elements. The characteristic roughness sizes k are selected so that 
they are a fixed and small fraction of the half-channel height h. The roughness Reynolds 
number k + is increased by increasing the overall friction Reynolds number Re T , which 
is the typical way roughness is encountered in practice (experimentally, and in practice, 
the roughness Reynolds number is usually increased by increasing the flow speed). 

In the Set-1 simulations (modelled roughness), the equations of motion, (12.ip . are 
numerically solved between two no-slip, impermeable walls at z = 0, 2 h and the e ffect o f 
roughness is represented by the parametric-forcing model of lBusse fe Sandham ( 2012 1. 
whereby a forcing term is added to the right-hand side of (12.11 a) of the form, 


— aF(z,k)u\u\5n, where F(z,k) = 


1, z < k or 
0, otherwise 


2h — k < z 


( 2 . 2 ) 


In principle, the effect of any roughness geometry, which includes both pressure and 
viscous drag, can be formally written as a body force on the right-hand side of the 
Navier-Stokes equation, but for the present purposes, the form (12.21) is meant to rep¬ 
resent a hypothetical homogeneous roughness. The roughness forcing is only active in 
the streamwise direction, opposes the mean flow and always dissipates kinetic energy. 
The roughness factor, a , is t hought to scale w i th the roughness density , that is, the 
frontal area per unit volume ( Nikora et al\ 2007 : Busse fe Sandhainl 2012h . measured in 
inverse-length (area per unit volume) units. Presently, k = h/40 and a = l/(40fc) for all 
simulations in Set 1 (refer to tabic Q] for details). A cosine mapping is used to stretch the 
grid in the wall-normal direction. 
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Figure 2. Grid for the Set-2 Re r ~ 180 minimal-span simulation. 



= 141 


In the Set-2 simulations, the sinusoidal roughness previously studied by IChan et~cH. 
( 20i nil is explicitly represented by a roughness-conforming grid composed of hexahedral 
cells (figure O • The grid is stretched in the wall-normal direction, with the stretching 
ratio held below 1.05. Grid information for all cases in Set 2 is detailed in table \B The 
roughness surface described by the sinusoidal function, 

z w = k cos(2tt x/X x ) cos(2n y/X y ), (2.3) 

locates the no-slip, impermeable wall relative to the reference plane at z = 0. Under this 
definition, we consider the characteristic roughness height k to be the semi-amplitude of 
this sinusoid. For the present sinusoidal roughness, the peak-to-valley roughness height k t 
(= 2 k), the root-mean-square roughness height k rms (= k/2) or other physical roughness 
heights could instead be used. Ultimately, the present method seeks to determine k s /k 
directly and the specific choice of the physical roughness height is therefore irrelevant, 
as all one requires in practice is the ratio between the equivalent sand-grain roughness 
and some readily measurable geometric property of the rough surface. For the present 
sinusoidal roughness, k s /k = 2k s /k rms = k s /kt/2. The roughness surface at the top wall 
is the mirror image (across the centreline) of the roughness surface at the bottom wall. 
Presently, k = h /18 and X x = X y « 7.1 k for all simulations in Set 2 (refer to tableQ]for 
details). 
























































Figure 3. Mean velocity profiles of simulated turbulent channel flow at (a) Re T « 590 and (6) 
Re T « 4000: solid, smooth; dashed, modelled roughness; grey, minimal span; black, full span. 
The vertical dashed line marks the top of the roughness-forcing region. The inset shows that the 
velocity shift stays constant for both minimal- and full-span channels above the roughness-forcing 
region. The various grey profiles in (6) correspond to the various minimal spans, Ly « 116, 232 
and 463. The inset in ( b ) shows that Ly « 116 is too narrow to accurately capture the roughness 
function for roughness heights above k + = 100. The full-span Re T « 4000 data in ( b ) is from 
iBernardini et a/I(120141) . 


3. Results and discussion 

3.1. Modelled roughness 

In figure [31 a), we show the mean velocity profile from four simulations with modelled 
roughness (Set-1 simulations) at Re T ~ 590, corresponding to the smooth full-span, 
the rough full-span, the smoot h minim a l-span and the rough minimal-span simulations. 
Consistent with the findings of iHwana ( 2013li . figure [3j a) shows that the minimal span 
of L+ ss 116 is sufficient for capturing the mean smooth-wall velocity profile given by 
the full-span simulations up to the height z+ ss 40. We will use the notation z+ to 
refer to the point above which the minimal profile departs from the full-span profile. 
The present simulations show that this is also the case for the profile above modelled 
rough walls (given by the dashed lines in figure Oa). Above the roughness height of 
k + = Re T /{h/k ) ~ 590/40 « 15, the inset of figure [3|a) shows that Af7 + for the full- 
span channels is reproduced by the minimal-span channels, even though U + itself above 
z+ ss 40 for the minimal-span channels fails to capture U + of the full-span channels. 
We will use the notation A U + to refer to the nominally constant value of the downward 
velocity shift at matched Re T , Uf —1/+, evaluated above the roughness height. The simi¬ 
larity of the mean relative velocity is, of course, a manifestation of Townsend’s outer-layer 
similarity: the downwards shift in velocity is purely determined by the near-wall flow over 
the roughness, which sets the wall drag. This shift will remain the same irrespective of 
the state of the outer profile, which is vastly different for minimal- or full-span channels. 
These results demonsrate the efficacy of the present method. It is important to empha¬ 
sise that the aforementioned method to evaluate A U + requires matched Re T , which is 
simple to achieve computationally since the outer length scale and the driving pressure 
gradient (and therefore U T ) can be fixed between smooth- and rough-wall simulations. 
Experimentally, where Re T between smooth and rough walls are seldom well-matched, 
A U + is typically evaluated from a shift in the log region (provided it exists). 

When the top of the roughness-forcing region is higher than z+, a wider minimal 
channel is required. A systematic study with various channel spans at Re T ss 4000, as 
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Figure 4. The Hama roughness function at fixed h/k = 40 for the modelled roughness, where 
k s /k ss 1. 6 , sim ulated in channels. The superimpo sed uniform sand-grain roughness data of 
iNikuradsel ( T93.il ) are extracted from I Jimenez; ( 2004 ). 


listed in tabled studies this effect. These results are shown in figurc[3](&). For the smooth 
walls, the data suggest that the minimal channel is able to reproduce the full channel 
profiles up to zf ~ 0.4L+. That is, the minimal profile from L+ « 116 reproduces the 
full profiles up to ss 0.4 (116) ~ 46, the minimal profile from L y ss 232 reproduces the 
full profiles up to z+ « 0.4 (232) rs 93, the minimal profile from L+ « 463 reproduces 
the full profiles up to z+ « 0.4 (463) « 185, and the pattern presumably goes on as 
long as z^~ ~ 0.4 L+ remains in the log layer z + < 0.15 Re T , where the size of eddies 
is thought to sca l e wit h distance from the wall. This behaviour is first reported by 
iFlores fc .Timenej ( 2010 ). who suggested Z+ ~ 0.3 T+. We refer to the region z < z c as 
being ‘unconfined’ since the turbulence in this region is faithfully captured, showing no 
signs of being constrained by the minimal span. A difficulty arises when the top of the 
roughness rises above the unconfmed region (when k > z c ). For the narrowest channel, 
we observe that the roughness function is over-predicted, giving AU + « 11 instead of 
the A U + « 9 obtained for the two wider (yet still minimal) channels. This sensitivity 
leads to an additional constraint for the minimum channel span, both of which must be 
satisfied 


L y > fc/0.4 and L+ > 100 (3.1) 

Physically, this means that the channel must be wide enough to fully immerse the rough¬ 
ness elements ( z c > k ) in natural, unconfined wall turbulence (which only occurs when 
Ly > 100). Although the constraint that L y > fc/0.4 is developed from the present mod¬ 
elled homogeneous roughness, where fc is the modelled roughness height, we expect that, 
in general, L y > O(fc)/0.4 for any physical roughness height fc such as k rms or k a ■ In 
practice, a convergence study with wider (still minimal) channels can be performed (cf. 
figure [3]&) to manage this effect. A conservative approach is to take L y > k t / 0.4, where 
fc t is the maximum peak-to-valley roughness height. The third constraint, requiring that 
L y > X y (where X y is some characteristic spanwisc length scale of the roughness texture, 
or repeating unit of the roughness) is not active for this case of modelled roughness, 
where the homogeneous forcing suggests X y —» 0. 

A sweep in roughness Reynolds numbers, from fc + ss 5 to 100, corresponding to Re T « 
180 to 4000 (tablelTJ is needed in order to fully characterise the roughness transition from 
the hydraulically smooth to the fully rough regime. This in turn enables determination 
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Figure 5. Mean velocity profiles of turbulent channel flow: (a) Re T ps 180, smooth and explicit 
roughness (h/k = 18) in minimal- and full-span channels; (6) Re T ps 180, 360, 540, 1080, smooth 
and sinusoidal roughness ( h/k = 18) in minimal-span channels. The vertical dashed line marks 
the top of the roughness element. The inset shows that the velocity shift stays the same for both 
minimal- and full-span channels above the roughness elements. 


of the equivalent sand-grain roughness k s /k = C , where C is a constant that best scales 
the highest roughness data onto the fully rough asymptote. Through the use of minimal- 
span channels, even the simulations at Re T ~ 4000, are feasible. The Hama roughness 
function for this sweep is presented in figureQ] The data from the full-span simulations at 
Re T ~ 590 and 950 are shown by the filled circles. It is noted that the roughness function 
profile obtained from the minimal-span simulations matches very well with the full-span 
simulations. Fitting to the fully rough regime, shown by the solid curve on figure [4} the 
present data show that k s /k ps 1.6 for a = l/(40fc) and so we have characterised the 
present modelled roughness from the transitionally rough regime, kf ps 7, to the fully 
rough regime, kf ps 160. The obtained value for k s could now, in theory, be used to predic t 
full-scale performance under this roughness condition ( Granville] Il958l : Schultz! 20071. 
Comp aring the minimal data with the uniform sand-grain roughness data of INikur adse 
( 19331) shown in light grey markers, reveals that the parametric-forcing model with the 
box-profile shape function closely mimics the rou ghne ss tra nsition of uniform-sand-grain 
roughness (as previously noted bv lBusse fe Sandham 2012 ). 


3.2. Sinusoidal roughness 

Figure[5](a) is a plot of the mean velocity profile from both the full and the minimal L+ ps 
141 channel with explicitly represented sinusoidal roughness elements (Set-2 simulations) 
at Re T ps 180. Consistent with the simulations with modelled roughness, the mean profiles 
from the minimal channel collapse onto the mean profiles from the full channel below 
about zff ps 0.4 L+ ps 56. For z + > 56, the minimal profiles diverge from the full channel, 
but the inset in[5])a) shows that the roughness function, AU + , remains the same for both 
full and minimal channels. In other words, the minimal channel is sufficient for predicting 
the roughness function. The method works in this case because the top of the roughness 
elements at z + = k + ps 10 is well within the unconfmcd region z + ps 0.4 ps 56, and 
the roughness wavelength (ps 71) is smaller than L+. An issue that often arises is the 
appropriate wall-normal location to evaluate the roughness function. In high-Reynolds- 
number flows, this location is frequently taken to be somewhere in the log layer. Presently, 
we observe that Uf~ — Uf~ is relatively insensitive to this location as long as A U + is 
evaluated above the roughness height. The insets in figures [3] and [5] demonstrate that 
this observation is consistent for all our results, for both explicit and modelled roughness. 
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th e maximum block age ratio, R/k = 6.7 5. The superim posed uniform sand-grain roughness data 
of iNikuradsel (119331 1 are extracted from I Jimene d (12004 ). 


Figure[5| b) shows the mean velocity profiles for the minimal channel at a fixed blockage 
ratio h/k = 18, where k + is increased from 10 to 60 by increasing the friction Reynolds 
number from Re T s=s 180 to 1080. For the lowest k + « 10, the span of the minimal 
channel is chosen to be L+ m 141 in order to properly capture the near-wall cycle 
and accommodate two periods of the sinusoidal roughness. In the cases with the higher 
roughness Reynolds numbers, k + ~ 20, 30 and 60, the minimal channel only captures one 
period of sinusoidal roughness, that is, X v = L y . The corresponding predicted roughness 
functions A U + are shown in the inset of figure [5]) &). 

In order to assess the efficacy of the present method for an explicitly represented 
sinusoidal roughness, we plot the Hama roughness function versus the equivalent sand- 
grain roughness kf in figure [6] Fitting the data in the fully rough regime, we are able 
to characterise this particular sinusoidal roughness with X x = X y s=s 7.1 k to find that 
k s /k sa 4.1. The minimal-channel predictions, shown by the open ma rkers, collapse ont o 
previously studied full pipe-flow simulations of the same roughness bv ldran et all ( 2015 ). 
albeit at various pipe blockage ratio R/k, as shown by the filled circles on figure [6] This 
collapse further validates the ability of the minimal channel simulations to recover the 
same estimate for k s as given by full-size simulations. 

It is worth pointing out that recent DNSs have also b een simulated in the fully rough 
regime ( Busse fe Sandham 20121 : lYuan fe Piomelli 2014 ). An important distinction here 
is that these previous cases used domain sizes that are deemed sufficiently large to capture 
the full range of turbulent motions. In contrast, the present minimal simulations have 
a severely confined outer flow with much lower intensities and a radically altered wake 
profile. The potential saving in computational cost for the minimal simulations could 
be redeployed to explicitly represent the roughness (see the present Set-2 simulations), 
without the need for roughness models. 


4. Conclusions 

We have presented a novel, fast and direct method for characterising the hydraulic resis¬ 
tance of any given surface roughness. The way in which a particular roughness transitions 
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from the hydraulically smooth to the fully rough regime is, to the first approximation, 
described by how the roughness geometry interacts with the near-wall flow. The method 
presented herein shows that this interaction, so far as the mean drag is concerned, is ac¬ 
curately captured using minimal-span channels. There are other roughness effects, such 
as the change in turbulent kinetic energy, and it remains to be seen the how these ef¬ 
fects are represented by the minimal-channel technique. We have validated the present 
method with a specific sinusoidal roughness and showed that the method is able to fully 
characterise the hydraulic (drag) behaviour of roughness not only in the transitionally 
rough regime, but also in the fully rough regime, yielding an estimate of k s which closely 
matches that given from full-span simulations. 

The dynamic drag characterisation of a rough surface is encapsulated in the Hama 
roughness function, which we show can be accurately determined using minimal-span 
channels. The fact that A U + is plotted against k + (or kf ) in the literature and not 
against Re T acknowledges that A U + , to a large extent, depends only on the roughness- 
affected near-wall flow. The minimal-span channel is a method for simulating this near¬ 
wall flow of thickness 0 (k) without resolving the outer scale h 3> k, thereby breaking the 
curse of the (outer) Reynolds number. The savings in computational cost for the present 
method are possible because capturing only the near-wall flow requi res far less gr id 


points than capturing the full flow. For example, the full-span DNS of IBernardini et al 
(12014 at Re T ~ 4000 requires N x x N y x N z = 8192 x 4096 x 1024 = 3.4 x 10 10 
grid points. In contrast, the minimal-span (L+ ss 463) DNS at the same Re T (table 
[TJ, which is sufficient to reach the fully rough regime, requires only N x x N y x N z = 
384 x 96 x 1024 = 3.8 x 10' grid points, or a three-orders-of-magnitude reduction in 
number of grid points. The actual wall-clock time will, of course, depend on the code 
and parallelisation, as well as the machine. As an indicative comparison using the same 
code and machine, the full-span ( L y /h = tt) Re T « 950 simulation (table [l]) requires 
2600 CPU hours/(/i/C/ T ), while the minimal-span (L+ « 112) Re T « 950 simulation 
(table [T]) requires only 42 CPU hours/(/ i/£/ t ). These figures amount to a two-orders-of- 
magnitude reduction in CPU hours. 

The present method can be used to characterise the drag characteristics of many sur¬ 
faces very quickly. Such a powerful tool enables the researcher to now focus on how the 
geometry of roughness affects the turbulent flow instead of focussing on setting up expen¬ 
sive and time-consuming experiments, physical or numerical. Indeed, the present method 
even enables researchers to reassess the physics underlying the success and failures of 
previously proposed geometrical factors. Recall that the minimal-span method retains 
the same benefit from a full-span direct numerical simulation in that the friction velocity 
formed from both viscous and pressure contributions is directly measured. 

In some ways, the present idea is not entirely new. In the large-eddy simulation subgrid- 
scale modelling methodology, the geometry-dependent large eddies are directly simu¬ 
lated whilst the subgrid small eddies are assumed to be universal and are modelled or 
understood through the energy-cascade phenomenology of Richardson, Taylor and Kol¬ 
mogorov. Presently, this reasoning is reversed for the roughness problem, in which the 
geometry-dependent small eddies are directly simulated whilst the universal ‘supergrid’ 
large eddies are assumed to be universal and are modelled or understood through the 
outer-laye r si milari ty hyp othesis of Townsend, a kind of small-eddy simulation, a term 
coined bv ljimenezl ( 20031 ). 

We have also provided guidelines on selecting the minimal channel. To obtain an ac¬ 
curate prediction of A U + , the roughness element must be submerged within unconfined 
near-wall flow, which yields the first two constraints: L y > k/0.4 and L+ > 100 re¬ 
spectively. Potentially, if k were very large, the first of these constraints could lead to 
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increasingly large ‘minimal’ channels. However, in practice it seems unlikely that these 
simulations would need to be conducted for 100 in order to establish the fully 

rough regime, and hence L y defined by this constraint should remain manageable. More 
problematic is the constraint that L y > X y . For heterogeneous or sparse rough surfaces, 
the spanwise length scale or repeating unit X y could potentially be very large in terms of 
viscous scaling, presenting an obvious Achilles’ heel for this proposed method, although 
this would be equally true for full span simulations. Even for surfaces which have a small 
X y , yet where the roughness is randomly arranged (such as a sanded or painted surface) 
the minimal span would need to be sufficiently large to contain a statistically representa¬ 
tive sample of the rough surface. Even so this X y would in most cases be significantly less 
that the irh box width often required of full-span simulations. Within these constraints, 
this method could potentially be extremely useful for all rough surfaces that exhibit ho¬ 
mogeneity over a relatively small scale. Examples could include surface finishes obtained 
from machining, painting or spray coating. One could also potentially use minimal chan¬ 
nels scaled in this way to investigate drag reducing rough surfaces such as riblcts. Indeed 
other passive or active surfaces where the primary effect of the surface is an alteration 
of the near-wall structure, or near-wall slip, could also be investigated using minimal 
channels (e.g. super-hydrophobic surfaces, compliant walls, porosity, acoustic liners). 

We can also point to certain possible improvements or extensions of this technique 
in the future. One issue with the minimal channels is that the centreline velocity be¬ 
comes very high (relative to U T ), which adds computational expense through required 
reductions in the time step to satisfy the Courant-Friedrichs-Lewy (CFL) condition. A 
possible solution here would be to add a body-forcing term at the channel centreline to 
appropriately manage the magnitude of this velocity. An alternative is to use a simple 
eddy-viscosity term that is active sufficiently far above the roughness elements to ac¬ 
count for absent eddies that are wider than the minimal span, which would reduce the 
centreline velocity. It is also possib le that the number of grid points could be further 
halved by using an open channel (Scotti 2006h . From a different approach, it could be 
the case that wider minimal channels may be more efficient in terms of obtaining sta¬ 
tistically converged statistics for given CPU hours (due to reductions in the ‘burstiness’ 
associated with minimal channels where L+ is close to 100). Finally, and perhaps most 
promising, we are keen to explore the capability of a single temporal simulation, starting 
from the highest Re T , and reducing the bulk velocity to map A U + versus k + (and hence 
obtain fc+) in a single numerical experimental sweep. 

This research was undertaken on the NCI National Facility in Canberra, Australia, 
which is supported by the Australian Commonwealth Government and also on the Victo¬ 
rian Life Science Computational Institute (VLSCI). The authors would like to gratefully 
acknowledge the financial support of the Australian Research Council (ARC). 
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